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Abstract 

A theory of the probability distribution function (PDF) tails of the 
blob density in plasma edge turbulence is provided. A simplified model 
of the fast convective radial transport is used. The theoretically pre- 
dicted PDF tails corroborate earlier measurements of edge transport, 
further confirming the strongly non-Gaussian feature of edge trans- 
port. It is found that increasing the cross sectional spatial scale length 
(L x and L y ) of the blob results in larger transport whereas increasing 
the toroidal scale length (L z ) decreases the PDF. The results imply 
that the PDF decreases for larger blob speed v\>. 
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I Introduction 



It is well known that turbulent transport determines the confinement of plas- 
mas in magnetic fusion devices. These inherently nonlinear phenomena are 
rather complex and still not well understood since they vary from improved 
confinement regimes to very violent disruptions. One important observation 
from experiments is that overall improved confinement is found when the 
edge turbulence is suppressed. Edge plasma turbulence is crucial for wall 
erosion and plasma contamination through the transport of particles and 
heat to the vessel walls, and thus for the confinement in future reactors pQ- 
[2]. Interestingly, experimental measurements of edge turbulence has shown 
the highly intermittent nature [3]- [5]. Furthermore simulations of statistical 
properties of edge turbulence in 2D j6]- [8] and 3D [9]- [11] have revealed 
generic non-Gaussian probability distribution functions (PDF) of fluctuation 
levels. In particular, in the turbulence simulation a large fraction of large 
events (or blobs) that ballistically propagate in the radial direction have been 
observed [12]- [U], which are especially dangerous for confinement. 

The so-called blob is a coherent structure with a higher density than the 
surrounding plasma which is localized in a plane perpendicular to the mag- 
netic field B while extended along the field line. When a charge dependent 
drift such as those induced by curvature or centrifugal force is present, the 
blob becomes polarized as the effective sheath resistivity creates an electric 
field. The resulting E x B drift transports the blob to the outer wall. The 
natural outward convective transport of blobs in edge plasmas indicates that 
these coherent structures may play a crucial role in intermittency and the 
non-Gaussian statistics in edge plasmas [2]. 

Coherent structures such as blobs, streamers or vortices are often asso- 
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ciated with avalanche like events of large amplitude and can therefore be of 
great importance for transport dynamics. Although these events are rela- 
tively rare, they can carry more than 50% of the total fluxes [16]- [17]. Con- 
ventional methods to characterize transport have been limited to mean field 
theory, where the transport is described by one averaged coefficient. There 
are however at present a lot of evidence that transport often involves events 
of many different amplitudes or scales, some of which are intermittent and 
bursty in time. Since these intermittent and bursty events are highly non- 
linear phenomena contributing to the non-Gaussian structure of the PDF 
tails, they are poorly described by mean field theory. To characterize the 
intermittent turbulent transport a non-perturbative way is needed [18]- [22] . 

In this paper we present a non-perturbative analytical theory of the PDF 
tails of density blob formation in tokamak edge plasmas. By adopting a sim- 
ple nonlinear fluid theory of the blobs [23]- [25], we predict that the blob 
density PDF tails have the exponential dependency (P(n b ) ~ exp{— (nf}), 
where rib is the blob density and ( is a coefficient dependent on the blob 
properties. Note that this scaling is similar to what was found for zonal 
flow structure formation in ion-temperature-gradient mode turbulence [22] . 
Interestingly, this exponential scaling agrees rather well with previous ex- 
perimental results with reasonable values of the coefficient ( for parameter 
values typical of plasma blobs [2]. 

Furthermore, we have found that increasing the cross sectional spatial 
scale length (L x and L y ) of the blob results in larger transport whereas 
increasing the toroidal scale length (L z ) decreases the PDF. Interpreting the 
constants in the blob speed seems to indicate that the PDF decreases for 
larger Vb- Note however that the blob speed is not a fundamental parameter 
in our model, but rather a combination of other parameters. 
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The paper is organized as follows. In Sec. II the physical model of 
the blob density is presented together with preliminaries of the path-integral 
formulation for the PDF tails of structure formation. In Sec III the instanton 
solutions are calculated and in Sec IV the PDF tails of blob formation are 
estimated. We provide numerical results in Sec. V and a discussion of the 
results and conclusion in Sec. VI. 

II Non-perturbative calculation of structure 
formation PDF 

The derivation of the physical model for radial plasma transport closely fol- 
lows Ref. [23]. We assume that the Scrape-off-Layer (SOL) plasma temper- 
ature, T, is constant. The electrostatic potential (f> is constant along the 
magnetic field B and can be calculated from the equation for electric current 

vjl + V||ii| = 0, (1) 

with j± = c(B x VP)/B 2 , where P = nT, n is the plasma density and c is 
the speed of light. Performing an integration along the field line and using 
the boundary conditions j\\\target = erifC s ^ at the targets we find 



e0 pi 



dlVlnB- (B x Vn) (2) 



T 2n t B 

where we have assumed that |^| < 1; n t is the plasma density at the targets, 
c s = \JT /M is the sound speed, M is the ion mass, e is the electron charge, pi 
is the ion gyro-radius, and the coordinate / goes along the magnetic field line. 
For a plasma blob with density rif, with parallel length L z situated around 
the midplane, Eq. (2) gives 
e0 _ L z pi dn b 

T ~ 2Rn t dy' U 



where we have neglected the magnetic shear and used VlnB = e x /R; R 
is the major radius and x and y are the local coordinates along the radial 
and poloidal directions, respectively. By using Eq. (3) for the E x B drift 
velocity, we find the blob plasma continuity equation in the form 

dn b CsPjLz ( d d 1 dn b . d . d 1 dn b S\ _ 
dt + 2R {dx [Ub dy [ n t dy )l dy [Uh dx [ n t dy } J ~ 1 [ ) 

When n t = for a constant £, the separation of variables gives 

n b (t,x,y) = n Q (x,t)e- {y/L ^\ (5) 
reducing Eq. 4 to a ballistic equation for n 

^ + ^^^°^'^ = °' 

with 

/ pi\ 2 L z n b 

= [Tj Rn t - (7) 
Note that the separable solution does not set the radial scale of the blob. 
The forcing / is defined in Eq. 8. 

There has been suggestion from both simulations and experiments is that 
the blob is formed from the non-linear saturation of the linear instabilities 
at the plasma edge [2j. Note that in the formation zone an approximately 
equal amount of enhanced density blobs and holes are generated [2]. The 
effective gravity (polarization) causes these newly formed coherent structures 
to move, the blobs move outwards whereas the holes move inwards. The 
ballistic equation describing the dynamics is symmetric under the change of 
parameters (x — > —x, no — > —no) [2E]. Note also that the blob velocity v b 
in Eq. (7) changes sign under this transformation. A detailed mechanism 
for the source of blobs is outside the scope of the present paper. In the 
following, we thus simply assume that there is a stochastic forcing (e.g. due to 



instabilities) and investigate the likelihood [probability distribution function 
(PDF)] of blob formation triggered by this forcing. Due to the stochastic 
forcing, blobs become short-lived in time, as shall be seen later. 

In order to calculate the PDF tails of blob formation, we utilize the 
instanton method [27] • To this end, the PDF tail is expressed in terms of a 
path integral by utilizing the Gaussian statistics of the forcing / [27J. We 
assume the statistics of the forcing / to be Gaussian with a short correlation 
time modeled by the delta function as 

(f(x,t)f(x',t')) = 5(t-t')K(x-x'), (8) 

and (/) = 0. The delta correlation in time was chosen for the simplicity of 
the analysis. In the case of a finite correlation time the non-local integral 
equations in time are needed. Note that the forcing / was chosen to excite 
blobs; the source of the forcing is assumed to be the fluctuations due to 
instability. The spatial overlap between the forcing and the blob is critical 
for the generation of a blob. 

The probability distribution function of blob density n b can be defined as 

P{Z) = (5(n b -Z)) 

dX exp(iAZ)(exp(— iXn b )) 

d\exp(i\Z)h, (9) 



where 

I x = (exp(-iAn 6 )). (10) 

The angular brackets denote the average over the statistics of the forcing /. 
The integrand can then be rewritten in the form of a path-integral as 

Vn b Vn b e~ s \ (11) 



where 



• f J2 m- ( drib , c sPi L z ( d d 1 dn b d d 1 dn b \ 

-,y d xdtn b [- W + ^ F - (^K^(-^)i - a^(^) J 



) 



+ 



+ 




(12) 



Note that P(Z) represents the probability of blob density taking a value Z. 

Note that the PDF tails of blob density can be found by calculating the 
value of S\ at the saddle-point in the case A — > oo. This will be done in Sec. 



Ill Instanton (saddle-point) solutions 

We have now reformulated the problem of calculating the PDF to a path- 
integral in Eq. (9). Although the path integral cannot in general be calcu- 
lated exactly, an approximate value can be found in the limit A — > oo by 
using a saddle point method to compute PDF tail. Since a direct application 
of the saddle-point equations results in very complicated partial differential 
equations for n b and n b , we assume that the instanton saddle-point solution is 
a temporally localized blob. That is, we assume that a short lived non-linear 
blob solution exists to the system of Eq. (4) in the form of a ballistically 
traveling solution Eqs (5) and (6). The blob density instanton takes the 
form n b (x,y,t) = no(x,y,t)F(t) while the target density is assumed to be 
n t = C n o- Here n (x,y,t) = n (x — v^e^^^ denotes the spatial form of 
the coherent structure or blob and F(t) is a temporally localized amplitude, 
representing the creation process. 

The action S\ consists of three different parts; the blob model, the forcing 
and structure formation, respectively. The full action including the forcing 
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+ 



Here, 



N 
N 



and structure formation terms can then be expressed in terms of the time 
dependent function F and the conjugate variable N, 

+ -kqJ dtN 2 

iXN J dtF5(t). (13) 

J d 2 xn {x - v b t)e-^ L ^ 2 , (14) 
J d 2 xn b (t, x)n (x - v b t)e~ {y/Ly) \ (15) 

K = (16) 

and the radial scale-length (L x ) is defined as 

dn 1 , 

-dx- = ~T x n °- (17) 

Kq in Eq. (13) is the strength of the forcing function k(x — x'), which is 
approximated by Taylor expansion in x and x' for simplicity. Keeping only 
the zeroth order terms in the expansion gives us the separable integral in 
x and x' . The time dependent function iV is the mean value averaged over 
the blob, N is the conjugate variable acting as a mediator between the real 
variable (N) and the forcing (/) and K is a constant used to simplify the 
expressions. 

The saddle point equations for instantons (the equations of motion) are 
obtained by minimizing the effective action S\ with respect to the indepen- 
dent variables F and N: 

w = -i- s+ r/ +K wrt FR )- XNm = °- (19) 



The equation of motion for F is derived for t < using Eqs. (18)-(19) as, 
ldF 2 



2 dF 

where 



= 77^ + 6772^ + 8773^, (20) 



m = (f ) 2 , (21) 

V2 = JJJJ^ = v 7 ^' ( 22 ) 

" 3 = W (23) 

Note that the constants 771, 7/2 and 773 all have the dimension of l/[Time] 2 . 
We first perform an integration in F and then use separation of variables to 
find 

r dF 1 r dF 

J FVr7i+ 477 2 F + 477 3 F^ " F(F + g) 

fM^x^i) = /* = «. (24) 
We then solve for F to find 

= FT^ exp{ ^ t} ' (25) 

F(0 = ^L_^%_. (26) 
W 2 m l -H(t) 1 J 

Note that H(t) < 1 and in Eq. (26) we impose the boundary conditions 

F(t -> -00) -> and F(0) = F . The initial condition for F(0) = F is 

found by integrating Eq. (19) over the time interval (— e, e) 

i(N(e) - N(-e)) = XN, (27) 

where we use the fact that the conjugate field N disappears for positive time 
(t > 0), which gives a relationship between the conjugate field N and the 
large factor A 

N(-e) = iXN. (28) 



The property of the conjugate variables to vanish for (t > 0) can be inter- 
preted as a causality condition. Using Eq. (18) at t — 0, we obtain 

— i^F+ j^F + 1 ^KF^j = -k N = -i\k N. (29) 

In the large A limit this gives 

F = ±-^=JL 2 y L x ^ N^/\. (30) 
2\/K v 

The instanton solution is localized within a time interval proportional to 
1 / The time scale of the blob or instanton solution can be estimated as 
r b iob = Vi = = 10~ 4 s, by using values in Ref. [2j. Comparing this blob time 
scale with that of the ambient fluctuations (r tur b = 10" 5 s) we find that T tUT b < 
Tbiob- This suggests that blobs described by this model with lifetime longer 
than the turbulence may significantly contribute to intermittent phenomena 
in the tokamak edge. 



IV The PDF tails 

We will now compute A dependence of S\ in Eq (13) at the saddle point i.e. 
the saddle point action which will then determine the PDF 



S x = -i J dtN [F + j-F + T ,P 
+ iXN J dtF5(t) + j dtN 



^! dt { p+ T: F+ '< F2 ) 2+ ' XNF ° < 31 > 



— / dt 4F 2 + i\NF (32) 
2k J 

^Lp^ + iXNF, (33) 



Where, 



7 = W' (34) 
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By using the initial condition for F [Eq. (30)], we find 

S x = -ia\ 3/2 (35) 



2 y3/2 L l L ^ K 

3 V K ' 
The PDF is found from Eq. (8) by utilizing the saddle point method 

P(Z) = J d\e- iXZ - s > (37) 



= J d\e- tXZ - iaXV2 (38) 

Let /(A) = — i\Z — ia\ 3 / 2 and find a A such that / attains its maximum 
and compute that value. This gives A = (§f ) an d 

/(*>) = < 39 > 

thereby giving PDF 

P(Z) = e^ z \ (40) 

C = 27^ < 41 > 
According to the definitions, for all reasonable physical situations the pa- 
rameters involved in ( are positive definite. Note that the saddle point so- 
lution justifies our assumption that A — > oo corresponds to Z — > oo. Eq. 
(40) provides the probability of finding a blob density (Z normalized by N). 
Note that when the forcing vanishes (/c — > 0, a — > 0) the PDF tails vanish 
(P(Z) - 0). 



V Results 

We have presented first prediction of the PDF tail of blob formation. By us- 
ing simplified model for the fast convective radial transport, we have found 
an exponential PDF tails of the form ~ exp{— C^l} similar to what was found 
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in Ref. [22] for zonal flow formation in ion-temperature-gradient mode turbu- 
lence. In this section the parameter dependencies of ( will be studied in de- 
tail. The results will be compared with a Gaussian prediction ~ exp{— (n 2 ,}. 

In Figure 1 (color online) the PDF tails of blob formation as a function 
of blob density {n b ) is shown by using the parameters 7 = jjj^j = 0.33 (red 
line, dash-dotted), 7 = 0.66 (blue line, solid line) and 7 = 1.35 (black line, 
dashed line) with k = 6.0. The Gaussian distribution (green line, dotted 
line) is also shown for 7 = 0.66 and k,q = 6.0. It is clearly shown that the PDF 
tails from the theoretical prediction recapture the experimental results shown 
in Ref. [2], where the PDF tails can be approximately fitted as ~ exp{— C n b} 
with T = 2.5 — 4.0. A decrease in the parameter 7 in Eq. (34) increases the 
PDF tails. The predicted PDF tails deviates significantly from the Gaussian 
distribution. 

To elucidate the parameter dependence of the constant ( a plot of l/a 2 
[a is given by Eq. (36)] as a function of the parameter 7 with kq = 3.0 
(black line, dash-dotted line), k,q = 6.0 (blue line, solid line) and k,q = 12.0 
(red line, dashed line) is displayed in Figure 2 (color online). Note that the 
PDF tails decrease as the parameter 7 decreases. In physical terms this 
means that increasing the cross sectional spatial scale length (L x and L y ) 
of the blob results in larger transport whereas increasing the toroidal scale 
length (L z ) decreases the PDF. Interpreting the constants in the blob speed 
seems to indicate that the PDF decreases for larger Vb- However, as shown 
in Sec. II the blob speed is not an independent parameter in our model but 
rather a combination of other parameters [see Eq. (7)]. In order to fit the 
theoretically predicted PDFs to experiments an estimation of the constant 
l/a 2 is needed. To this end, we used experimental values in Ref. [2] and 
obtained an estimate for the constant l/a 2 as 1.5 ± 0.5. 
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Figure 1: (Color online). The blob density PDF tail as a function of n& 
normalized by N. The parameters are 7 = jjj^j = 0.33 (red line, dash- 
dotted), 7 = 0.66 (blue line, solid line) and 7 = 1.35 (black line, dashed line) 
with k = 6.0. the Gaussian distribution (green line, dotted line) is shown 
with 7 = 0.66 and k = 6.0. 




Figure 2: (Color online). I /a 2 as a function of 7 with kq = 3.0 (black line, 
dash-dotted line), k = 6.0 (blue line, solid line) and k q = 12.0 (red line, 
dashed line). 
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VI Discussion and conclusions 

In order to elucidate the highly intermittent turbulent transport in the edge 
and SOL we have considered a non-perturbative model of intermittent trans- 
port driven by blobs. By using a non-linear model of fast convective edge 
plasma transport derived in Ref. [23], we utilized the instanton calculus [27] 
and [18]- [22] to calculate the PDF tails of blob formation. The result- 
ing PDF tails have the exponential form ~ exp{— (nf}. Interestingly, this 
prediction agrees with a number of experimental results showing a highly 
non-Gaussian statistics of the transport at the edge [2]. Furthermore, the 
PDF tails from experiments show a generic exponential form ~ exp{— (n^} 
with T = 2.5 — 4.0, corroborating our prediction. Note that this is the first 
calculation of PDF tails of blob formation, which were shown to be strongly 
intermittent. Considerable transport can however be mediated by rare events 
of large amplitude assuming that density blobs cause radial transport. 

It is also of interest to study the momentum flux driven by blobs. To this 
end, we replace blob formation (A / d 2 xdtn b 5(t)) in Eq. (13) by the blob flux 
/ d 2 xdt[n b -§-(±-^)}5{t)) to obtain the following action 



Using Eq. (42) to find the instanton solutions gives us similar results as 
before, but with the initial condition Fq oc A. This initial condition then 
makes the scaling for the PDF tails of momentum flux (here denoted by Z) 
as P(Z) ~ exp {— (^zZ 3 / 2 }. This is similar to what was found for momentum 
flux in drift wave turbulence [18]- |21j. 




(42) 
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It is interesting to note that the exponential scalings of the predicted 
PDFs are similar to those found for zonal flow formation and momentum flux 
in ITG turbulence. The reason for this ubiquitous exponential PDFs with 
the same scaling is because the order of the highest non-linear interaction 
terms in the governing equations is the same, giving the same dependence of 
the large parameter A in the initial conditions (Fq = F(0)), and thus similar 
exponential scalings of the PDF tails. The non-linear term can be easily seen 
to be quadratic in blob density from Eq. 4. 

We have shown that the PDF tails of blob formation depends on the char- 
acteristic scale lengths (L x , L y ,L z ). Although, there is no direct dependency 
on the blob velocity which is not an independent parameter in our model, it 
is indicated that the PDF tail is inversely dependent on blob speed (due to 
decreasing (). Moreover the size of the blob is crucial for the PDF, with a 
larger poloidal (L y ) and larger radial extension (L x ) giving larger PDF tail. 

We note that a non-Gaussian scaling of the PDF (the exponent of n&) is 
found even when the forcing is Gaussian, although the exact exponent may 
depend on the temporal and possibly spatial correlation of the forcing (/). 
In the present paper, the forcing is chosen to be temporally delta correlated 
for simplicity. The source of the forcing is assumed to be the fluctuations. 
In general it would be desirable to identify a relation between the forcing 
and the turbulence amplitude. However, this is still an open problem and 
unfortunately outside the scope of the present paper. 

Finally, we find a good agreement between our predicted PDFs and the 
experimental results reported in Ref. [2J. In particular the exponential scal- 
ing is very similar to that found in experiments for reasonable parameter 
values. These findings strongly suggest that the experimental results are due 
to intermittent phenomena coming from edge turbulence. In a future publi- 
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cations we will consider the issues of blob generation and the asymmetrical 
PDFs found in many experimental and numerical work. 
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